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^ . ABSTRACT 

Context. With the new generation of high-resolution X-ray spectrometers the understanding of warm absorbers in Active Galactic Nuclei has 
improved considerably. However, the main question remains the distance and structure of the photoionised wind. 

Aims. We study the absorption and emission properties of the photoionised gas near one of the brightest and most variable AGN, the Seyfert 
galaxy NGC 405 1, in order to constrain the geometry, dynamics and ionisation structure of the outflow. 

Methods. We analyse two observations taken with the Low Energy Transmission Grating Spectrometer (LETGS) of Chandra. We study the 
spectra of both observations and investigate the spectral response to a sudden, long-lasting flux decrease of a factor of 5 that occurred during 
the second observation. 

O I Results. We confirm the preliminary detection of a highly ionised component with an outflow velocity of -4500 kms"', one of the highest 
velocity outflow components seen in a Seyfert 1 galaxy. The sudden drop in intensity by a factor of five during the second observation causes a 
drop in ionisation parameter of a similar magnitude in the strongest and main ionisation component (v = -610 km s"'), allowing us for the first 
time to determine the recombination time of this component and thereby its distance in a robust way. We find an upper limit to the distance of 
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, 10 m, ruling out an origin in the narrow emission line region. In addition, an emission component producing strong radiative recombination 
• • 1 continua of C vi and C v appears during the low state. This can be explained by emission from an ionised skin of the accretion disk at a distance 
_ ' of only ~ 4 X 10'^ m from the black hole. Finally, the spectra contain a broad relativistic O vm line with properties similar to what was found 

, before in this source with XMM-Newton; this line has disappeared during the low flux state, consistent with the disappearance of the inner part 
» I of the accretion disk during that low flux state. 

Ci Conclusions. Combining high-resolution spectroscopy with timing information, we have constrained the geometry of the emission and ab- 
sorption components in NGC 4051. 

Key words, galaxies: Seyfert - quasars: individual: NGC 4051 - galaxies: active - X-rays: galaxies 



1. Introduction 

The stupendous amount of energy emitted by an Active 
Galactic Nucleus is released by gas that flows towards a su- 
permassive black hole in the centre of a galaxy, presumably 
via an accretion disk. The average energy of the emitted ra- 
diation increases towards the centre of the accretion disk, and 
thus X-rays provide the best probe of the gas flow in the im- 
mediate surroundings of the black hole. From the presence of 
jets, as well as from blue-shifted ultraviolet absorption lines, 
we learn that in addition to the flow towards the black hole, 
there is also gas flowing away from it. With the availability of 
X-ray spectrographs on Chandra and XMM-Newton, we can 
study these flows close to the black hole, by studying lines in 
the X-ray spectra. The presence of multiple absorption line sys- 
tems, that differ in their level of ionisation or in their outflow 
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velocity, c onstrains the geome try of the ou tflow (for example 
NGC 3 783. lNetzer et alJll20()l : NGC 5548 ■ ISteenbrugge et alj 
il200.5h l. 

In this paper we describe the X-ray spectra of the Active 
Galactic Nucleus of NGC 4051, obtained with the Low- 
Energy Transmission Grating on board of the Chandra satel- 
lite. Emission from the nucleus of NGC 4051 was already 
found by Hubble (1926), and the galaxy is one of the 6 ob- 
served by ISevfertI li 19431) . and discussed in his seminal pa- 
per as one of the first 12 known Seyfert galaxies. Being a 
member of the Ursa Major cluster jTuUv etalJI 19961) it has a 
distance of 18.6 Mpc JTuIIv & Pierce "2000') and a redshift of 
700 kms ' llVerheiien& Sancisi 200n . NGC 4051 is a rela- 
tively bright optical object, at V ^ 13.5. In X-rays, the source 
was not detected with UHURU or Ariel- V, but was studied 
extensively with EXOSAT and ROSAT; its X-ray variabil- 
ity was recently analysed on the basis of RXTE and XMM- 



2 



M. Fenovcik et al.: High-resolution X-ray spectroscopy of NGC 4051 



■ 1 ■ 1 

A 


■ 1 ■ 1 ■ 
B 

I'M 

' Jan 2002 - 


1 . 1 . 1 . 1 . 



20 



40 60 
time [ks] 



80 



100 



Jul 2003 



.-)-■*-*-++++ 



20 



40 60 
time [l<s] 



80 



100 



Fig. 1. The light curve in the zeroth spectral order (top panels) and hardness ratio (bottom panels) for the observations of Jan 
2002 and July 2003. The hardness ratio is defined as the ratio of the 2-8 and 8 - 32 A counts. The top panel has a bin size of 
350 s and the bottom panel has a bin size of 2000 s. Dotted lines denote the different flux states during the observations. 



Newton data (^McHard v et"ai]l2004 . Its X-ray continuum has 
been described with a power law with photon index F =1.8-2.0 
jNandra & Pounds II1994I) . A soft excess with respect to this 
power l aw is not wel l described by a multi-temperature disk, 
and Ogl e et alJ ( l2004 concluded that it is due to broad emis- 
sion from relativistic O vii emission lines and to radiative re- 
combination. 

Two absorption line systems, at -2340 + 130 and -600 + 
130 kms were found in the X-ray spect rum with the 
High Energy Transmission Grating of Chandra JCoUinge et alJ 
I2OOII) . Multiple absorption systems were also found in the ul- 
traviolet with the Space Telescope Imaging Spectrograph and 
in the Far Ultraviolet with the Far Ultraviolet Spectroscopic 
Explorer; one of these may coincide with the -600+1 30 km s ' 
X-ray abs orption syste m, all the others have sma ller absolute 
velocities jColUnge et al.-2001.:iKaspi et al.l2004 . 

Inter estingly, a preUminary an alysis of our first LETGS 
spectrum lvan der Meer et al ] (l2003h did not show that high ve- 
locity component, but indicated the presence of an outflow 
component at even higher velocity, -4500 km s ' . For that rea- 
son we requested our second observation. This second obser- 
vation, taken when NGC 405 1 was in a much higher flux state, 
allowed us to investigate the ionised outflow both in terms of 
outflow velocity and ionisation structure. In this paper we re- 
port the analysis of both observations with the LETGS. 

2. Observations and data reduction 

In this paper we describe two observations of NGC 405 1 with 
the Low Energy Transmission Grating Spectrometer (LETGS, 
iBrinkman et alJ lEoOO.) ) of Chandra. The detector used was 
the High-Resolution Camera (HRC-S). The observations were 
done from December 31, 2001 17:40:41 to January 1, 2002 
19:51:46 and July 23, 2003 00:34:15 to July 24 03:24:43, and 
with 94.2 and 96.6 ks net exposure time, respectively. Below 
we refer to these observations as the Jan 2002 and Jul 2003 ob- 
servati on, respectively. The data reduction is described in de- 
tail bv lKaastra et alJ (l2002i) . Basically, after event selection we 



produce background subtracted spectra and light curves. As the 
HRC-S detector has limited energy resolution, higher spectral 
orders cannot be separated in the data but are taken into account 
properly during spectral fitting. 

3. Data analysis and results 

For all spectral fittin g we used SPEX code developed at SRON 
jKaastra et alJl9 9^). Details of the spectral model and the time 
intervals for which spectral fits are made are given later in this 
paper We find that each spectrum requires several components 
to be fitted adequately. To minimise the number of components, 
we add them one by one, until the added component does not 
improve the fit significantly. It will be noted that the contin- 
uum emission components and the absorption components de- 
scribe the global spectrum, whereas the narrow and broadened 
line emission features are limited to specific wavelengths. We 
therefore determined errors on the continuum emission and ab- 
sorption systems while keeping the emission line parameters 
fixed to the values of the best fit; and determined errors on 
the emission line features while keeping the continuum and 
absorption lines fixed to the values of the best fit. The fit pa- 
rameters are listed in Tabled All models were corrected for 
Galactic ab sorption with a column density of 1.31 x 10^"* m"^ 
(lElvis et ai...l98 9) and also for the small cosmological redshift 
of 700 kms"^. Finally, Table |2] gives the increase in;t'^ when 
omitting a given spectral component from the final model. All 
errors correspond to A^^ = 2 (84 % confidence). 

3.1. Flux level variations 

NGC 405 1 was variable in both our observations, as shown in 
Fig.^ where we plot the count rate obtained with the LETGS 
in zeroth spectral order During the Jan 2002 observation the 
average brightness of the source increased about halfway the 
observation to a somewhat higher level; during the Jul 2003 ob- 
servation the source was brighter yet during the first ~ 65 000 s 
and then dropped in about 3000 s to the level of the second half 
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Table 1. The parameters of the spectral components for the January 2002 whole observation (plus A and B part) and July 2003 
whole observation (plus C and D part). The time intervals A to D are defined in Fig.^ All wavelengths and velocities refer to 
the rest frame of NGC 405 1 . 



January 2002 



A part 



B part 



July 2003 



C part 



D part 



Continuum 

• power-law 

Flux (2-10 keV)' 
Photon index T 

• modified black-body 
Flux" 

kT (keV) 



0.96: 
1.85: 



0.04 
0.04 



1.2 ±0.2 
0.102 ±0.003 



0.83 ±0.12 
1.58_02, 



1.0 ±0.3 

103+"°" 



1.19 : 
2.01 : 



:0.05 
:0.05 



1.2 ±0.2 
0.103 ±0.005 



1.92: 
2.24: 

2.2: 
0.143 : 



0.07 
0.03 

0.3 
0.005 



2.57 : 
2.23 : 



0.06 

:0.02 



2.8 ±0.3 
0.149 ±0.003 



0.41 : 
2.22 : 

1.0: 
0.179: 



0.03 
0.13 

0.2 
0.020 



Absorption components 

• warm absorber P 
Nh 

• warm absorber 2^ 
Nh 

V(ari, 

• warm absorber 3^ 
Nh 

log^ 

V(ari, 
Vo,„ 



0.34 ± 0.05 
0.67 ± 0.06 

150+^*" 
-210 ±20 



0.7 ± 0.3 
1.80 ±0.25 

75 ±30 
-620 ± 60 

19+15 

3.2 ±0.1 
160 ± 60 
-4450 ± 110 



1 95+220 
2.37 ± 0.26 

-680^«° 



74+74 

3 4+"' 



0.6 ± 0.2 
1.40 ±0.18 

90 ±40 
-570 ± 80 



0.6 
2.24 
90 



0.3 
0.10 

+40 
20 

570 ± 60 



14 : 
3.0 d 



0.1 



160'- 



-11230 

-30 

-43 10+™ 



130+™ 
-4570 ± 120 



20± 10 

3 r" ' 

1 j + 100 

-467o"'»'' 



1000 



0.41 ± 0.04 
0.74 ± 0.06 
150 ±20 
-210 ±20 

1.9 ± 1.0 
2.44 ± 0.08 

40 ±25 
-610 ±40 

15+'° 
3.2 ±0.1 
15 ± 15 
-4600 ± 120 



3.5 ± 1.0 
1.54 ±0.10 

50 ±20 
-580 ± 40 



Broad emission lines 
• O vm Lya line* 



flux 


12.7 ± 1.8 


9.7 ± 1.8 


15.9 ±2.6 


6.9 ±0.7 


11^1 


1.0 ±0.4 


X 

FWHM 


[21] 


18.5+°? 
[21] 


[21] 


19.26_oio 
[2.5] 


\9.Q%\\* 
[3.2] 


18.77 ±0.16 
0.69 ± 0.32 


q 

i 

• O vn triplet* 


5.2 ± 0.6 

48+' 

^°-22 


5 2+'° 
49+13 


5.2 ± 0.6 

49+5 


1 9+0.2 

49+40 


2.0+°2 

46:{ 




flux 








3.8 ±0.7 


5 0+" 

-'•"-1.8 


1.4 ±0.4 


X 








22.07 ±0.14 


22.07 ±0.16 


21.72 ±0.02 


FWHM 








1.50 ±0.30 


1.61 ±0.40 


0.16 ±0.06 


• C VI Lya line* 














flux 








2.3 ±0.4 


2.4 ± 0.6 


1.8 ±0.5 


X 








33.68 ± 0.03 


33.70 ±0.05 


33.68 ± 0.04 


FWHM 








0.30 ± 0.08 


0.32 ±0.13 


0.25 ± 0.07 


Narrow emission features 














• Nil forbidden line* 














flux 


0.76 ±0.12 


0.9 ± 0.2 


0.8 ±0.2 


1.0 ±0.2 


0.9 ± 0.2 


1.4 ±0.2 


X 


22.089 ± 0.007 


22.089 ± 0.006 


22.087 ±0.013 


22.095 ± 0.005 


22.093 ± 0.006 


22.090 ±0.012 
















T 


3 2+2' 

■'•^-1.0 


7 Q+6.8 
^•^-14 


T 9+3.0 
-'•^-1.1 






4.9 ± 1.8 


EMcv 


0±20 


0± 100 


0±20 






180 ± 100 


EMcvi 


40 ±30 


20+3° 


70 ±50 






320 ± 120 


EMnvi 
EMnvi, 


10+-° 
0± 10 


10+40 

0± 10 


0±20 

7+16 

'-7 






30 ±40 
0±20 


EMovii 


7±7 


7+11 

'-7 


11 ± 11 






16± 15 


EMoviii 


0±5 


0±5 


0± 15 






10± 10 



' Absorption-corrected flux in 10"'* Wm"^ 

2 Absorption-corrected bolometric flux in lO"'"* Wm"^ 

3 The column density Nh is in units of 10^^ m"^. The ionisation parameter ^ is in units of 10"' Wm, the turbulent velocity 
velocity Vo,„ are in km s"' 
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V,,,,/, and the outflow 

The flux is in phm"'^ s"', the wavelength and the FWHM are in A. In the case of the relativistic profile, / is in degrees. FWHM values in [] 



were not fitted but computed from the best-fit relativistic line profile. 

^ the temperature is in eV and the emission measure is in units of lO'"* m"^ 



4 M. Fenovcik et al.: High-resolution X-ray spectroscopy of NGC 4051 

Table 2. Increase inx^ when a given spectral component is omitted from the model. 



omitted component 


number of 
parameters 


Ton ir\f\i 
J an zuuz 


A 
/\ 


D 


July zuuj 




u 


warm absorber 1 


4 








380.56 


413.94 




warm absorber 2 


4 


41.05 


7.06 


51.54 


102.08 


80.34 


55.58 


warm absorber 3 


4 


76.27 


34.12 


54.36 


59.89 


20.63 




VIII hya line 


4 or 3 


41.04 


25.76 


27.71 


54.34 


45.28 


7.51 


VII triplet 


3 








34.84 


45.28 


17.23 


C VI Lyo- line 


3 








22.79 


31.39 


12.38 


VII forbidden line 


2 


57.65 


36.94 


23.03 


43.74 


23.80 


29.37 


RRC 


7 


7.95 


3.66 


5.68 






30.92 


of full model 




1614.84 


1603.74 


1602.77 


2820.43 


2755.70 


2144.82 


d.o.f. 




1478 


1478 


1478 


2563 


2563 


1591 




Observed wavelength (A) 



Fig. 2. Flux spectrum and our best spectral model for each individual light curve part, which are marked as depicted on Fig. [T] 
Note the presence of strong RRCs of C vi at 25.43 A and of C v at 3 1 .80 A in spectrum D. The O vii forbidden line at 22. 10 A 
is visible in every spectrum and has the same strength in spectra A, B and C, but it is stronger in spectrum D. Note the different 
continuum level and shape in spectrum C compared to the other spectra. 



of our first observation. It then recovered gradually to the high 
level again. On top of these long-term variations we also detect 
vaiiability with smaller, but significant, amplitude on shorter 
time scales of a few thousand seconds. These rapid variations 
are smaller after the rapid flux decrease of the Jul 2003 ob- 



servation. In the analysis below we will discuss the spectra at 
different flux levels, defined by the lightcurve of Fig.^ A and 
B are the first and second half of the Jan 2002 observation, re- 
spectively, and C and D are the first and second interval of the 
Jul 2003 observation. Expressed in units of 10"'"^ Wm"^, the 
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0.5-2.0 keV fluxes of these spectra are 0.44, 1.11, 3.41 and 
0.66 for A, B, C and D, respectively. 

In Fig. ^we also show the hardness ratio, defined as the 
ratio of the count rates in the 2-8 A and 8-32 A band passes. 
The hardness ratio of the Jan 2002 observation, HR - 0.215(7) 
is higher than that of the Jul 2003 observation, HR = 0.160(3), 
where the number between brackets is the 1-sigma error in the 
last digit. The hardness variation between different flux levels 
is small, with HR = 0.239(14), 0.203(8),0.161(3), 0.134(21) 
for intervals A, B, C and D, respectively. 

3.2. Continuum 

In both observations the continuum cannot be fitted with only 
a simple power law. The soft excess can be modeled by 
adding a black body modi fied by coherent Compton scatter- 
ing (Ka astra & Barrl[l989l) . The parameters of the power law 
component of the Jan 2002 and Jul 2003 observations reflect 
the results found from the lightcurve and hardness ratio, which 
show that for the Jan 2002 observation the normalisation is 
three times smaller than for the Jul 2003 observation, while 
the slope is flatter, i.e. the spectrum is harder in Jan 2002 (see 
TablelD- 

In contrast the normalisation of the modified black body is 
the same within the error in both observations. The black body 
temperature is higher in the Jul 2003 observation (kT - 0.143 
keV, vs. 0.102 keV in the Jan 2002 observation). 

While looking at the different flux levels within each obser- 
vation, we find that the modified black body component does 
not change between intervals A and B, whereas the powerlaw is 
harder at the lower flux of A than at the higher flux of B. In con- 
trast, the change in flux level between C and D is reflected in 
the normalisations of the power law, and of the modified black 
body, whereas the slope of the power law is the same in C and 
D. 

3.3. The ionised absorber 

Thanks to the excellent spectral resolution of the LETG instru- 
ment the presence of a warm absorber in the spectrum is not 
controversial. We fitted the warm absorber with model XABS 
in SPEX. In this model, the ionic column densities are not inde- 
pendent quantities, but are linked via the ionisation parameter 
^ - L/nr^, where L is the source luminosity, n the hydrogen 
density and r is the distance from the ionising source. In addi- 
tion to ^, the fit parameters are the hydrogen column density of 
the absorber A^//, the outflow velocity v'o„, and the turbulent ve- 
locity of the medium Vmrb- The advantage of the XABS model 
is that all relevant ions are taken into account including also the 
ones with the weakest absorption features. 

3.3.1 . The warm absorber in the January 2002 
observation 

During the January 2002 observation the source is in a low flux 
state and consequently the S/N ratio is lower This has a nega- 
tive effect on the detection of the fainter absorption lines which 



cannot be detected. Only a few strong absorption lines can be 
used for the determination of the warm absorber properties. 

With XABS we find that two warm absorber components 
are needed. The component with the higher velocity (-4450 
kms ') has a higher ionisation parameter and a higher col- 
umn density than the component with the lower outflow veloc- 
ity (-620 kms"'; see TableQ. The absorption system with the 
higher outflow velocity improves the fit by Ax^ - 76 for 4 ad- 
ditional degrees of freedom. This is the fastest outflow velocity 
yet found in a Seyfert 1 galaxy. 

We fit the spectra A and B starting with the warm absorber 
components found in the total spectrum of Jan 2002. The sepa- 
rate spectra for A and B then are found to have within the error 
bars the same parameters for the warm absorber components 
as the total Jan 2002 spectrum. Due to the large uncertainty 
of the parameters for spectrum A, we are not very sensitive 
to detect changes; thus while the factor of 3 luminosity differ- 
ence between A and B could in principle give rise to a similar 
change of the warm absorber properties, we cannot detect such 
a change. 

3.3.2. The warm absorber in the July 2003 observation 

To fit the Jul 2003 observation we need three absorption sys- 
tems. Two of these are virtually identical to the components 
fitted in the Jan 2002 observation; the third component has a 
lower outflow velocity (-420 km/s), and also a lower ionisa- 
tion parameter and lower column density than the other two 
(see TableQ. 

In spectrum D we detected only one absorption system, at 
an outflow velocity of about -790 kms"'. The absorption col- 
umn is higher for spectrum D than the absorption column for 
the warm absorber component with the same velocity detected 
in interval C and in the total Jul 2003 observation; and the ion- 
isation parameter is lower. 

3.4. Broad emission lines 

Besides absorption lines also strong emission features are 
present in every spectrum. We first describe the broad emission 
lines. 

The presence of broad emission features is illustrated in 
Fig.|3]which shows excess flux above the power law and mod- 
ified black body continuum for the Jul 2003 observation. Clear 
excess emission in a band with a width of at least an A is 
present near 19 A and 22 A. This is the region of the ls-2p 
transitions in O viii and O vii, respectively. In addition, there is 
excess emission near 34 A, the region of the ls-2p transition in 
C VI. The Jan 2002 observation also shows a significant excess 
but only near the O viii transition. 

Close to the black hole the formation of an emission line 
is influenced by the strong gravity field and therefore the final 
emission line has a relativistic line profile. We modeled this 
relativistic profile as a narrow emission line convolved with the 
relativistic disk line profile of lLaor 1 (ll99ll) . The inner (ri) and 
outer (r2) disk radius were fixed in our models to the default 
values (ri = l.234GM/c^ and r2 = 400GM/c^). The disk inch- 
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nation / and the disk emissivity slope q (emissivity proportional 
to R '' with R the distance to the black hole) were allowed to 
vary. In the case when the emission line is created further away 
from the black hole this profile shows no or only a small devi- 
ation from the classical gaussian profile. 




Observed wovelength (A) 

Fig. 3. Data and continuum model (power law and modified 
black body spectrum) in the O vm Lyor region for the C part 
(upper panel) and the D part (lower panel). Here the absorption 
and emission line features are not yet taken into account. Note 
the broad excess near 19 and 22 A. 



In order to model the broad emission lines and to search for 
possible time variability in particular parts of the observation 
we proceeded analogously as in the case of the absorption fea- 
tures. First we determined the best model for the emission lines 
for the full observations and after establishing that emission 
model we used it as the first approximation for the emission 
model in particular parts. 

The parameters that we found for the broad lines are shown 
in TableC] 

3.5. Narrow emission features 

In both full observations we detect a narrow Ovii forbidden 
line. The Ovii forbidden line at 22.10 A is also visible in 
each individual spectrum and has the same strength in spec- 
tra A, B and C, but it is stronger in spectrum D. The line 
shows only a weak blueshift: the average outflow velocity is 
-110 + 50 kms-i. 

The main difference between spectrum D and the other 
spectra is the presence of strong Radiative Recombination 
Continua (RRCs) from C vi (25.30 A) and C v (31.62 A) in the 
spectrum, see Fig. Hand Tabled The temperature of the recom- 
bining photoionised gas is about 5 eV. In addition the RRCs 
show a redshift with respect to the rest frame of the galaxy 
of ~ 1000 kms We do not find significant detections for 
RRCs of Nvi (22.46 A), Nvii (18.59 A), Ovii (16.77 A) and 
O VIII (14.23 A). In all other spectra we do not detect significant 
RRCs, so the numbers given in Table ^ should be interpreted 
as upper limits. 



4. Discussion 

4.1. Tine sudden decline of the source flux in 2003 

In this section we focus upon the sudden decline of NGC 405 1 
in the July 2003 observation, at the boundary of part C and part 
D (Fig.[l). Within 3000 s, the flux dropped by a factor of ~ 5 
and remained low for more than 20 000 s. Large flux variations 
are common in this object. In fact, around 29 ks after the start 
of the Jul 2003 observation, also a drop in intensity of a factor 
of 5 occurred, but this was after a rise by a factor of 2.5 during 
the previous 2000 s, and as Fig.^shows, the recovery to higher 
flux level was much faster in that case. 

The uniqueness of the transition between part C and D is 
the rather abrupt change preceded and followed by long periods 
of relative quiescence. This makes it possible for us to see the 
response of the continuum spectrum, the broad line emission 
and the warm absorber to this sudden flux drop. 

The main change between C/D can be characterised as fol- 
lows (see also Table[0: 

1. The modified black body component becomes somewhat 
hotter but with three times less bolometric flux. If we as- 
sume that this component represents the direct emission 
from the accretion disk, this can be explained if the inner 
part of the disk has disappeared and therefore does not radi- 
ate anymore. The only modest change in temperature then 
may imply that the range of disk temperatures in the disap- 
peared part of the disk is not too different from that of the 
remaining part. 

2. The drop in power law flux is stronger, a factor of ~6, while 
the photon index does not change. The disappearance of 
the inner part of the disk leads to less seed photons for the 
assumed inverse Compton scattering corona, not only be- 
cause of the lower disk flux (~3 time less disk flux) but also 
due to the relatively lower fraction of disk photons from the 
outer, remaining parts of the disk that are captured by a hot 
corona that is close to the black hole. The fact that the pho- 
ton index hardly changes implies probably that apparently 
the optical depth and temperature of the corona remains the 
same 

3. The relativistic O vm line that can be seen in spectrum C is 
not visible in spectrum D. This also agrees with the scenario 
of sudden disappearance of the inner accretion disk. 

We therefore suggest that the sudden decline between part 
C and part D is caused by the disappearance of the inner ac- 
cretion disk. Such phase transitions are not uncommon in other 
galactic and extragalactic accretion disk systems. Whether the 
disappearing part of the disk is blown away through a strong 
wind, through jets or is being accreted onto the black hole is 
unclear. 

The precise fraction of the disk that disappeared is a lit- 
tle hard to estimate and strongly model dependent; we only 
measure that about 2/3 of the bolometric modified black body 
flux disappears. Around a maximally rotating disk, 2/3 of the 
the emitted disk flux comes from within ~ IGMjc^, but due to 
general relativistic effects (gravitational redshift, cone of avoid- 
ance causing most photons to be re-captured by the disk or 
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black hole) this constitutes much less than 2/3 of the observed 
disk flux. On the other hand, for a non-rotating black hole, with 
a smaller cone of avoidance, 2/3 of the flux originates from 
within ~ 2QGM/c^. Given this, we estimate as an order of mag- 
nitude that the disk within ~ IQGM/ c^ has disappeared. W ith a 
black hole mass of M = 3 x 10^ Mq i McHardv et alJ2004 ). the 
light crossing time for the diameter of the disappeared region 
is 30 s. We observe that the fading takes a full 3000 s, hence 
the average disappearance speed is only 3000 km s"' . 

4.2. The ionised absorber 

4.2.1 . The high ionisation component 2 

In our highest quality spectrum, part C of the July 2003 ob- 
servation, we detect three absorption systems (Table Two 
of these are at a relatively low velocity (components 1 and 2) 
while component 3 has a high velocity. 

The X-ray component 2, at an outflow velocity of - 
610 kms"', has a column density of ~ 2 x 10^^ m"^. Apart 
from several L-shell lines from iron ions, its dominant imprint 
on the absorption spectrum are the Lya transitions of hydro- 
genic oxygen, nitrogen and carbon. Column density, outflow 
velocity as well as average 2-10 keV flux of spectrum C are 
very similar to the values found by Pounds et al. ( 2004) in their 
analysis of the RGS spectrum taken in May 2001. Therefore 
we identify our component 2 with their dominant component. 
However, in our case the ionisation pa rameter 4^ is a facto r 
of 10 larger than the value obtained bv Poun ds et alJ (|2004'). 
Various processes may contribute to this. First, while the 2- 
10 keV flux in both these observations are the same, our pho- 
ton index is ~ 0.4 larger, leading to a much softer spectrum. 
For spectrum C the total ionising flux may be a factor of 4 
higher, leading to a much higher ionisation parameter for gas at 
the same distance and density. A further dilTerence is that dur- 
ing our observation the relative variability during observation C 
was relatively small (Fig.^, while during the RGS observation 
th e flux varied wil dly (see for instance the light curve shown 
bv iMcHardv et aU (.2004.)) . Ther efore, the ionisation parame- 
ter foundbylPound^^O (I2OO4I) might be an average between 
high and low states, if the warm absorber responds fast enough 
to the continuum variations. An indication for at least some 
spread of t he ionisation parame ter due to variability in the RGS 
data is that P ounds et al . ( 2004) also find a high ionisation com- 
ponent (log^ - 2.7) at the same outflow velocity. In our more 
"quiet" spectrum C, we do not find evidence for such a highly 
ionised component at the outflow velocity of -600 km s"' . 

Interestingly, when the continuum flux drops by a factor 
of five in ~ 3000 s, defining the boundary between spectrum 
C and D, the ionisation parameter of component 2 drops by a 
similar factor (8 ± 3). Because of the fortunate circumstance 
that the flux variations during intervals C and D were relatively 
small, we therefore conclude that the warm absorber responds 
rapidly to the flux decrease. 

How fast does the absorber respond? To that aim, we have 
plotted the transmission of component 2 for parts C and D 
(Fig. 15 central panel). The most striking features is a broad- 
band transmission decrease of 10-20 %. This by itself is 
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Fig. 4. Comparison between the model transmission of individ- 
ual warm absorber components for the C and D states. From 
top to bottom: components 1, 2 and 3. Dark curves: transmis- 
sion for the high state spectrum C based on the parameters from 
Table[n Light colour curves: for component 2, the transmission 
during the low state D; for components 1 and 3, the transmis- 
sion for the high state C but with ^ 5 times lower, therefore 
showing the expected transmission if each component would 
respond very fast to the factor of 5 flux decrease between C 
and D. The dotted lines indicate the 15-17 and 17-19 A bands 
used in Fig.|5l 
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Fig. 5. Upper panel: light curve in the 20^0 A band of the July 
2003 observation; bin size: 2000 s. Lower panel: fluxes in the 
15-17 A band (data points with error bars) and 17-19 A band 
(histograms; error bars similar to the other band); these narrow- 
band fluxes have been normalised to the average flux per A in 
the 20^0 A band shown in the upper panel. Therefore if there 
would be no spectral changes, the normalised 15-17 and 17- 
19 A fluxes should have remained constant. The dotted lines 
indicate the boundary of the low state D. 



not straightforward to measure in spectra with relatively poor 
statistics extracted for short time intervals, because a slight 
change in continuum parameters could mimic the transmission 
change. Also, the signal to noise ratio in individual lines is too 
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low. However, Fig. 0] shows that there is a sharp decrease in 
the 15-17 A band, caused by a deepening of the Ovn edge 
and the development of the Fe-M Unresloved Transition Array 
(UTA) complex between 15-16 A. Therefore we have made 
light curves in the 15-17 A band (Fig.|5}- A comparison with 
the 17-19 A light curve shows that both bands follow each 
other quite well during the high state C, but in the low state 
D the 15-17 A flux is on average 25 % lower than the flux 
above the O vii edge (17-19 A band), in agreement with Fig.|3] 

It is clear from Fig.|5]that the absorber responds definitely 
within 10 000 s to the continuum change; it may be much faster 
(down to the 1000 s scale or even faster), but the quality of our 
data does not allow us to prove that. 

We have estimated recombination time scales for arbitrary 
ions / from the general equation: 



dA^ 
"dT 



-NJi(L) - NiRiiT) + A?,-i//-i(L) + NmRm(T) 



«e (1) 



where A^, is the density of ion /, rif. the electron den- 
sity, Ri(T) the temperature-dependent recombination coeffi- 
cient from ion / to ion / - 1, and /,(L) is the luminosity- 
dependent ionisation coefficient from ion / to ion ; + 1 . For our 
order of magnitude estimate is is sufficient to neglect here any 
change in the shape of the ionising spectrum but to reckon only 
with the drop in luminosity L. We start with the high-state equi- 
librium (log^ = 2.44) and calculate the ion concentrations. We 
then decrease L by a factor of 5 and evaluate Q and determine 
the characteristic time scale r,- defined by d(lnA^,)/df = 1/t,-. 
In the high state (C), component 2 has an equilibrium temper- 
ature (for a stationnary state) of 20 eV. In the low state (D), 
component 2 has an equilibrium temperature of 5 eV. We there- 
fore determine t, for temperatures between 5 and 20 eV. We 
find that O vii has a characteristic time scale t, of lQ*/ni2 s for 
r = 20 eV, and 4 x 10Vni2 s for T = 5 eV, where nn is the 
hydrogen density in units of lO'^ m""*. For other ions, in partic- 
ular Fe XIII - Fe xvii, the main contributors to the Fe-M complex 
between 15-16 A in spectrum D, the time scales are 0.2 to 2.5 
times the time scale for O vii. 

We conclude that from our observed upper limit of the re- 
combination time scale of f^c ^ 10"* s, it follows that the den- 
sity n S 10'^ m"^. Using the definition of the ionisation param- 
eter (^ - Ljnr^), we then obtain an upper limit r ^ 10'^ m for 
the distance of the warm absorber of component 2 to the central 
source. From its measured column density of ~ 2 x 10^^ mT^, 
we find then an upper limit to the thickness d of the absorber 
oft/ ^ 2x 10" m. 

The distance upper limit of 10'^ m implies that we can rule 
out an origin in the narrow emission line region for component 
2. An origin a t the in ner edge of a dusty torus (^ 3 x 10'^ m, 
cf. lElvis et al. ' ( 2005)) cannot be fully ruled out but is unlikely 
given the fact that our recombination time is probably faster 
than 10"* s). Therefore, an origin in the outer parts of an accre- 
tion disk seems more appropriate. The relative thickness of the 
wind d/r ^ 2x 10""* is extremely small; this suggests that we 
have a very narrow structure, or if the outflow covers a larger 
region, it must have a small volum e filling factor. 

We note that in a recent paper ( lElvis et al]l2005h a recom- 
bination time of only 3000 s is proposed based on broad-band 



fits of time-resolved, low resolution EPIC spectra. That detec- 
tion lacks sufficient proof because it is based upon a correla- 
tion study of derived ionisation parameter versus count rate; at 
low spectral resolution, however, it is not evident ho w to disen- 
tangle possible rapidly changing relativistic lines jOele et al] 
ll2004 found such lines in the same XMM-Newton dataset) 
from a rapidly changing warm absorber In our case, we use 
high-resolution spectra and have a much cleaner transition from 
a high to a low state. 



4.2.2. Low ionisation component 1 

At only a slightly smaller outflow velocity (-200 km s"') com- 
pared to component 2, we find our second, lowly ionised com- 
ponent in spectrum C (Tabled component 1). This component 
is most easily seen through the deep ls-2p absorption lines of 
O vn, O VI and O v. The outflow velocity of both components 1 
and 2 coincide \ yith th e broad troughs seen in the FUSE Ovi 
Une lKaspj et'ai] ( l2004 and in the range of UV absorption c om- 
ponen ts seen in several lines observed by STIS^ollinge et alJ 
(1200 1! ). Component 1 has five times smaller column density 
than component 2, and is in good agreement with the low ion- 
isation tail seen in the column density versus ionisation plot of 
iQgleet alJ ( E004i) . 

We cannot see component 1 in the low state spectrum D, 
after the drop of a factor of 5 in intensity. The five times lower 
flux level combined with the shorter exposure time makes de- 
tection of any absorption component in spectrum D harder as 
compared to spectrum C. Further, even if the gas of this com- 
ponent responded similarly to the continuum decrease as com- 
ponent 2, it would be hard to find: Fig.l^shows that the contin- 
uum opacity in both situations is low for component 1 , and also 
there is no deep Fe-M UTA through that could be recognised 
after substantial binning of the data. We therefore cannot make 
any statement whether component 1 remained present, or that 
it got a lower ionisation. 



4.2.3. Higli-velocity ionisation component 3 

A high velocity component at an outflow velocity of 
-2340 kms ' was first found bv lCollinge et alJ I2OOI ); we do 
not confirm that detection but instead find a much higher out- 
flow v elocity component at v = -4500 kms ', as reported 
fist bv lvan der Meer et alJ (l2003h . This component appears to 
be present in both our spectra of 2002 and 2003, with al- 
most the same column density and high ionisation parameter 
(log^ - 3.2). In general, this component has a low opacity. Its 
predicted response due to the flux decrease between part C/D 
in the July 2003 observation is shown in Fig. |3 If the warm 
absorber in this component would be of sufficient low density 
(we estimate n ^ 10'"* m"^), it would respond rapidly to the 
flux decrease. That would result in a significant decrease of the 
transmission of about 20 % below the O viii edge around 14 A. 
We do not see clear signatures of this in our spectra, but this 
may be hidden due to the lower opacity than component 2, and 
perhaps also a part of the enhanced opacity (if present) could 
be mimicked in our spectral fits by a slightly different slope 
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for the power law component in spectrum D. In the absence of 
clearer evidence, the distance of component 3 remains unclear. 

4.3. Broad emission lines 

While in the observation of July 2003 we detect with > 3cr 
significance three broad emission features (O viii Lya, the O vii 
triplet and C vi Lya) in the observation of January 2002 only 
one broad emission feature is present (the Oviii Lya line). 
Adding a broad C vi Lya line or the O vii triplet does not im- 
prove the fit significantly. 

There is a distinct difference between the line profile of the 
broad O viii Lya line in both observations. In the observation of 
January 2002, the line is extremely broad (FWHM 21 A, with a 
sharp blue cut-off near 17 A and an extended red tail. This line 
profile is caused by the steep emissivity law (q - 5.2) com- 
bined with the small inner radius of the disk (corresponding to a 
rapidly rotating Kerr black hole). The emissivity slope q = 5.2, 
the inner disk radius r, and the disk inclination / = 50deg 
are in excellent agreement with the values given by Osle et al. 
( |2004 . These authors observed NGC 4051 with the RGS in- 
strument of XMM-Newton in May 2001, half a year before 
our observation. There is only a small difference in equivalent 
width: we find 4.7 A while Ogle et al. find 2.6 A. We conclude 
that in general the parameters of the relativistic line as found by 
us in the Jan 2002 spectrum are rather similar to the parameters 
found by Ogle et al. 

However, Ogle et al. modeled the soft excess in NGC 405 1 
with the help of relativistic emission from the entire O viii emis- 
sion Unes series and RJiCs. We tried to model the soft excess 
in the RGS (2-37A ) spectral band with the same model used by 
lOgleet al.l(l2004 (i.e., ignoring the modified black body com- 
ponent). This trial was not successfully neither for the January 
2002 observation nor for the July 2003 observation (see Table 
|3}- We conclude that in both LETGS observations the soft ex- 
cess requires also the presence of the modified black body 
emission, rather than the higher Lyman transitions of O viii. 

During the observation of July 2003 the profile seems to be 
produced under different conditions. We tested the possibility 
that also in this observation the O viii Lyo' profile is modified 
by relativistic effects. The inclusion of a relativistic broadened 
profile in the model is not statistically more significant than a 
simple gaussian profile. However, a gaussian model predicts a 
velocity broadening of the line of the order of 30000 kms ' 
and this velocity broadening may hide some relativistic smear- 
ing. The best fit for a relativistic profile is achieved with an 
emissivity slope q - 2, with a FWHM of 2.5 A. Thus the line 
profile is an order of magnitude narrower than in the obser- 
vation of January 2002. A similar behaviour of the emissiv- 
ity slope q, as a function of the central source flux, was found 
also for anothe r narro w line Seyfert 1 galaxy (MCG -6-30-15, 
llwasawa et al.l (Il996l) ). but in that case for the Fe-K line. 
Considering the same scenario as for MCG -6-30-15, the red 
tail during January 2002 observation of NGC 4051 suggests 
that during that low state the line should be produced very close 
to the black hole where gravitational effects are significant. The 
blue horn is gravitationally shifted into the red wing and almost 



Table 3. Comparison of the soft excess models for the January 
2002 and July 2003 observation. 



Observation 


model 




dof 


January 2002 


relativistic lines model" 


698 


561 




our best model* 


684 


565 


July 2003 


relativistic lines model" 


1640 


1148 




our best model* 


1326 


1152 



" the relativistic lines model consists of a power law plus th e entire se- 
ries of O vm lines broadened by the relativistic profile as in lOgle et alJ 
j2004h 

* our best model consists of a power law, modified black body and a 
single relativistic O viii Lya line. 

all of the line emission goes to the red wing. Conversely, during 
a high state most of the emission comes from a larger distance. 
Thus the blue peak comes most ly from the blue side of the disk 
and the red wing is depressed (" ( llwasawa et alJl996l) V 

The spectrum in July 2003 is dominated by the high state 
C. Surprisingly, after the sudden drop in intensity that occurs 
at the end of spectrum C, an extremely broad line as seen dur- 
ing the lower flux states A and B does not re-appear. Instead, 
there remains only a weak, relatively narrow broad emission 
line (FWHM 0.69 A or 1 1 000 km s '), which can be attributed 
to the outer accretion disk or the broad line region. 

The fact that the O viii Lya line may be explained with a 
relativistic profile while the O vii triplet has a gaussian profile 
does not to need be contradictory. Most of the O vm Lya line 
can be formed by the gas which is closer to the black hole and 
the gas producing the O vn line may be placed further away. 
However an important caveat is that O vii triplet shape is more 
difficult to determine as it is a blend. 

Relativistic line profiles in the soft X-ray band have been 
discovered first bv Branduardi-R avmont et al.. (,2001) in some 
other narrow lines Seyfert 1 galaxies (MCG -6-30-15 and 
Mrk 766). In these sources, the disk line parameters derived 
from these lines and the Fe-K line are in good agreement. The 
XMM-Newton observation of NGC 405 1 that showed a strong 
relativistic O viii line (Os\& et al. 2004) did not show a signif- 
icant relativistic Fe-K line (Pounds et^al, 2004). Unfortunately 
the LETGS is not sensitive enough in the the Fe-K band and 
thus we cannot directly compare the disk geometry as derived 
from iron and oxygen lines. For the physi cal implications of the 
relativistic lines in NGC 4051 we refer to lOele et 311(12004 . 

4.4. Narrow features 

During our observations we detect various narrow emission 
features in the spectra. The forbidden O vii line and RRCs of 
C V and C vi were the most dominant. 

The O VII forbidden line is narrow and does not respond to 
the flux change (Table 0. This suggests that this line is pro- 
duced by gas away from the central source, for example in the 
narrow line region. 

There is no evidence for any variation of the centroid of the 
O VII forbidden fine. Combining all observations, the centroid 
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of 22.093 + 0.004 A as measured in the restframe of NGC 405 1 
corresponds to a net outflow velocity of - 110 + 50 km s ' . 

The strongest RRCs are detected during a short time inter- 
val in the July 2003 observation (Fig. The shape and cen- 
troid of individual RRCs is not only affected by the temperature 
of the recombining gas and the emission measure of the rele- 
vant ion, but also by the motion (velocity broadening) of the 
recombining gas. In a previous part of our paper (section l33l 
we determined the temperature and emission measure assum- 
ing zero velocity broadening. How can the velocity broadening 
affect the determination of the temperature and the emission 
measure? We have investigated this by including in our RRC 
spectral model a velocity broadening component. On the basis 
of the goodness of fit (x^) we exclude that the recombining gas 
is moving faster than 3 000 km s"'. The temperature and emis- 
sion measure are within the error bars the same for models with 
a velocity broadening of 3 000 km s"' and for models without 
velocity broadening. 
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Fig. 6. Data and continuum model in the C vi RRC region for 
the first half (upper panel) and the second half (lower panel) of 
interval D (see Fig.Q of our observation of NGC 4051. 



If we interpret the velocity broadening as a typical 
Keplerian velocity at the location where the RRC is formed, 
we can derive some properties of th e recombining gas. For a 
black hole mass of M = 3 x IO^Mq ( McHardv et al.l2004 ) and 
the maximum velocity broadening of Av - 3000 km s"', we 
typically get a distance of the recombining gas of R ~ 4x 10'^ 
m or ~ 4500 Schwarzschild radii. 

From the fitted temperature of the RRC during part D (4.9 + 
1.8 eV) we derive an ionisation parameter log^ = 1.6 ± 0.3. 
With the ionising luminosity of the source (L = 2.8 x 10^^ W), 
our distance derived above and the definition of ^ we obtain 
a typical density of 4 x 10'^ m"^. It is now also possible to 
estimate the volume of the RRC emitting source. Using the ob- 
served emission measure of the C vi RRC of 3.2 x 10^^ m"^, a 
cosmic carbon abundance and a fraction of 50 % of the Carbon 
in the form of Cvii (appropriate for log^ = 1.6), we find 
an emitting volume of V = 1.2 x lO-'^ m-'. Equating this to 
V = QR^AR with AR the characteristic thickness of the emitter 



in the radial direction and Q. the solid angle sustained by the 
emitter, we find 

n— = 0.002 (2) 
R ^ ' 

Another constraint follows from the column density A'h 
through the emitting region. This can be written as A^h - 
nbJi - nR{Q.QQ2/Q.) or substituting more numbers A^h - 
2.5 x 10^^(47r/f2) m^-. Using that estimate, we can rule out 
an origin of the RRC emission in the warm absorber: the ioni- 
sation parameter log ^ ^ 1 .6 would suggest an association with 
warm absorber component 2, but the column density of that 
component is at least 10 times smaller. Instead, the high (col- 
umn) density may suggest an origin at the thin, ionised skin of 
the accretion disk at a few thousand Schwarzschild radii from 
the black hole. 

We conclude with some consistency checks. First, the light 
crossing time through the emitting region is of order 10^ s, con- 
sistent with the duration of phase D. Note that we see the RRC 
during the full interval D (Fig.|6j. Further, for log^ - 1.6 we 
expect only a small fraction of gas recombining into O vm, but 
instead we expect a strong O vii RRC at 16.77 A, but this fea- 
ture is not clearly observed. However in this spectral region 
absorption from ions such as Fe ix may mask the presence of 
O vn RRC emission. 

5. Conclusions 

In this study we have investigated the spectral properties of 
NGC 4051, observed by Chandra-LETGS on two occasions 
for a total of 180ks. 

The time averaged spectrum of both the Jan 2002 and the 
Jul 2003 can be fitted by a modifyed black body and a power 
law, absorbed by ionised gas. 

Two of the ionisation components (2 and 3) are well vis- 
ibible both in Jan 2002 and in the high state (part C) of Jul 
2003. These gas components are consistent with being stable 
over a long time scale (~ 19 months) both in ionisation level 
and outflow velocity. 

In particular, we report the detection of the highest outflow 
velocity-gas ever observed (v ~ -4800km s"'), for compo- 
nent 3. The lower ionisation phase of the gas (component 1) is 
well detected only in the state C. This lack of detection is most 
probably due to the lower statistics affecting the other time seg- 
ments. 

In the last part of the Jul 2003 observation, (part D), only 
the second component is detected, but it shows a significant 
variation in the ionisation parameter which linearly responded 
to the continuum flux variation on a time scale >3000s. From 
this we estimated a lower limit for the gas density (« ^ 
lO'^m"-') and as a consequence a distance for the absorber 
r $ 10'^ m and a thickness of the gas layer c/ ^ 2 x 10" m. 
The emission spectrum is rich in narrow and broad emission 
lines: 

The narrow lines (e.g. the O vii forbidden line) are not vari- 
able in time and therefore consistent with being produced in 
regions very distant from the black hole (e.g. the NLR). 
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The RRCs of C v and C vi showed instead a very rapid vari- 
ability, significantly arising above the continuum after the flux 
drop in spectrum D. The estimated distance (r ~ 4 x 10'^ m) of 
the recombining gas shows that these RRCs are produced close 
to (or in) the accretion disc. 

We also find that the column density of the gas produc- 
ing the RRCs is incompatible with that found for the warm 
absorber, suggesting a different location for the absorbing and 
emitting media. 

Broad lines of O vm, O vn and C vi are also evident in the 
spectrum. The O vn and C vi lines show a simple Gaussian pro- 
file with a FWHM ranging from 0.3 to 1.5 A, suggesting an 
origin in the BLR. 

On the contrary, the O viii Lya line is best modeled by a 
relativistically broadened profile. The tentatively detected vari- 
abiUty on a months time scale is according to expectations for 
a line produced this close to a black hole. In accordance with 
previous studies, we find that the observed anticorrelation be- 
tween the emissivity slope and the central source flux can be 
explained in terms of a relativistic accretion disc around a spin- 
ning black hole. 
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